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Abstract 

The dominant electroweak two-loop corrections to the precision observables Mw 
and sin 2 9 e g are calculated in the MSSM. They are obtained by evaluating the two- 
loop Yukawa contributions of 0(a 2 ), O(cttab), 0(a 2 ) to the quantity Ap. The result, 
involving the contributions from Standard Model fermions, sfermions, Higgs bosons 
and higgsinos, is derived in the gauge-less limit for arbitrary values of the lightest 
CP-even Higgs boson mass. A thorough discussion of the parameter relations enforced 
by supersymmetry is given, and two different renormalization schemes are applied. 
Compared to the previously known result for the quark-loop contribution we find a shift 
of up to +8 MeV in Mw and —4 x 10 -5 in sin 2 9 c fi . Detailed numerical estimates of the 
remaining uncertainties of My/ and sin 2 9 c r from unknown higher-order contributions 
are obtained for different values of the supersymmetric mass scale. 
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1 Introduction 



Electroweak precision observables (EWPO) like the masses of the W and Z bosons, Mw,z, or 
the effective leptonic weak mixing angle, sin 2 9 e s, are highly sensitive probes of the quantum 
structure of the electroweak interactions. The Standard Model (SM) and any extension or 
alternative predicts certain relations between these observables that can be tested against 
the corresponding experimental values. The experimental resolution is better than the per- 
mille level [1], and thus the measurements can be sensitive to even two-loop effects. Hence 
the EWPO are very powerful for discriminating between different models of electroweak 
interactions and for deriving indirect constraints on unknown parameters such as the masses 
of the SM Higgs boson or supersymmetric particles, see Ref. [2] for a review. A recent analysis 
[3] considered the EWPO, combined with the anomalous magnetic moment of the muon and 
the decay b — > 57, in the constrained MSSM with cold dark matter bounds. Already at the 
current level of experimental and theoretical precisions the sensitivity of the EWPO to the 
scale of supersymmetry allows to infer interesting information, pointing towards a relatively 
low scale of supersymmetric particles. 

The current experimental accuracies, obtained at LEP, SLC and the Tevatron, for M w , 
M z and sin 2 # cff are 5M W = 34 MeV (0.04%), 5M Z = 2.1 MeV (0.002%)and 5sm 2 9 cS = 
16 x 10~ 5 (0.07%) [1]. At the GigaZ option of a linear e + e~ collider, a precision of SMw = 
7 MeV [4,5] and 5sm 2 6 cS = 1.3 x 10~ 5 [5,6] can be achieved. 

The evaluation of the SM theory predictions has reached a high level of sophistication. 
The one-loop results for the M^-Mz mass relation, parametrized by the quantity Ar, and 
for sin 2 9 e s are completely known. For Ar also the full two-loop contributions are known [7,8], 
while for sin 2 9 e s the calculation of all two-loop contributions involving a closed fermion loop 
has recently been completed [9,10]. Leading universal corrections to the EWPO in the SM 
and extensions of it enter via a quantity called Ap. It parametrizes the leading universal 
corrections from vector boson self energies induced by the mass splitting between fields in 
an isospin doublet [11]. Within the SM various two- and three-loop corrections have been 
obtained [12-19]. A contribution to Ap induces the following shifts in M w and sin 2 eff (with 
\-sl^cl = M^ w lMl): 

AM W « ^^yA/,, Asm 2 eff «-^-\Ap. (1) 

In the minimal supersymmetric extension of the SM (MSSM) [20], the theoretical evalua- 
tion of the EWPO is not as advanced as in the SM. In order to fully exploit the experimental 
precision for testing the MSSM and deriving constraints on the supersymmetric parameters, 
it is desirable to reduce the theoretical uncertainty of the MSSM predictions to the same 
level as the SM uncertainties. So far, the one-loop contributions to Ar and sin 2 9 cS have 
been evaluated completely [21,22]. In the case of non-minimal flavor violation the leading 
one-loop contributions are known [23]. At the two-loop level, the leading 0(aa s ) corrections 
to Ap [24] and the gluonic two-loop corrections to Ar [2,25] are known. 

In the present paper we calculate the two-loop MSSM-corrections to the EWPO that 
enter via Ap at 0(a 2 ), 0(a t ab), 0(al). These are the leading two-loop contributions in- 
volving the top and bottom Yukawa couplings and come from three classes of diagrams with 
quark/squark loop and additional Higgs or Higgsino exchange (sample diagrams for the three 
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Figure 1: Sample diagrams for the three classes of contributions to Ap considered in this 
paper: (q) quark loop with Higgs exchange, (q) squark loop with Higgs exchange, (H) 
quark/squark loop with Higgsino exchange. denotes Higgs and Goldstone boson exchange. 



classes are shown in Fig. 1). These contributions are of particular interest, since they involve 
corrections proportional to m\ and bottom loop corrections enhanced by tan (3, the ratio of 
the vacuum expectation values of the two Higgs doublets of the MSSM. Partial results have 
already been presented in Ref . [26] . 

As a first step, in Ref. [27] the O{of), 0(a t ab), 0(al) corrections were calculated in the 
limit where the scalar quarks are heavy, corresponding to taking into account quark/Higgs 
diagrams (class (q)) only. While this class of corrections turned out to be well approximated 
by the SM contribution (setting the Higgs-boson mass of the SM to the value of the CP-even 
Higgs-boson mass of the MSSM), a potentially larger effect can be expected from diagrams 
with squarks and higgsinos, classes (g), (H) in Fig. 1, which do not possess a SM counterpart. 
In the related case of similar two-loop corrections to — the squark contributions turned 
out to be much more important than the quark contributions [28]. 

For the two-loop Yukawa corrections in the SM it turned out that the dependence on the 
Higgs-boson mass is numerically important. While the Higgs-boson mass is a free parameter 
in the SM, the masses of the CP-even Higgs bosons of the MSSM are given in terms of the 
other parameters of the model. In the "gauge-less limit" that has to be applied in order 
to extract the leading two-loop Yukawa corrections, the mass of the lighter CP-even Higgs 
boson, M/j, formally has to be put to zero. In Ref. [27] it was observed for the calculation of 
the diagrams of class (q) that can be set to its true value instead of zero in a consistent 
way. In the present paper we provide a detailed discussion of the gauge- less limit, yielding 
an explanation of this observation. We will analyze the Higgs-mass dependence also for the 
other classes of diagrams in Fig. I. 

We analyze the numerical effects of the new corrections for various scenarios in the 
unconstrained MSSM and for SPS benchmark scenarios [29]. We study two different renor- 
malization schemes and investigate the possible effects of unknown higher-order corrections 
for Myy and sin 2 e g . 

The outline of the paper is as follows. First we review the theoretical status of the EWPO, 
focusing on the role of Ap and the gauge-less limit (Sect. 2), then we describe the structure 
of the calculation (Sect. 3). In Sect. 4 the renormalization is discussed, taking into account 
the implications of the gauge- less limit. The numerical analysis is performed in Sect. 5. We 
conclude with Sect. 6. 
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2 Existing corrections to Ap in the MSSM 



The quantity Ap, 

Sz(0)_E^ 
P ~ Ml Ml ' {2) 

parametrizes the leading universal corrections to the electroweak precision observables in- 
duced by the mass splitting between fields in an isospin doublet [11]. Hz,w($) denote the 
transverse parts of the unrenormalized Z and W boson self-energies at zero momentum 
transfer, respectively. The corresponding shifts in the predictions of Myj and sin 2 9 e g have 
been given in eq. (1). 

A detailed overview of the existing corrections to the EWPO within the SM and the 
MSSM can be found in Ref. [2]. Here we briefly review the existing corrections to Ap in the 
MSSM and the corresponding SM contributions. Since Ap as defined in eq. (2) is not a phys- 
ical observable but contains only a part of the quantum corrections to the EWPO, some care 
has to be taken in its evaluation. Only contributions that form a complete subset of leading 
corrections to the EWPO can be consistently incorporated into Ap. While the fermion loop 
contributions to Ap form a meaningful subset of the full one-loop corrections to the EWPO 
(see below), a naive evaluation of eq. (2) for the bosonic one- loop contributions of the SM 
would result in an expression that is neither UV-finite nor gauge-parameter independent (see 
the discussion in Refs. [2,30]). A consistent inclusion of the bosonic corrections therefore 
makes it necessary to go beyond the Ap approximation. We will introduce the so-called 
"gauge-less limit" in Sect. 2.2 and show that at 0(a 2 ), 0(a t ab), C(ct 2 ) the contributions to 
Ap coincide with the complete contributions to the EWPO. 



2.1 One- loop and higher-order QCD corrections to Ap 

The dominant one-loop contribution to Ap within the SM arises from the top/bottom dou- 
blet [11]. It is given by 



3G 

A Pl-loop = - 

c 

where 



Ap? M ioo P = ^^^oK,m 2 ), (3) 



2 x y x 

F (x,y)=x + y log-. (4) 

x-y y 

F Q has the properties Fo(m 2 ,m|) = -^(m 2 ,, m 2 ), F (m 2 ,m 2 ) = 0, F (m 2 ,0) = m 2 . One 
therefore obtains F (m 2 ,m 2 ) ~ m 2 , giving rise to the well-known quadratic dependence of 
the one-loop corrections to the EWPO on the top-quark mass. 

Within the MSSM the dominant correction from supersymmetric particles at the one- 
loop level arises from the scalar top and bottom contribution to eq. (2). For m h ^ it is 
given by 

A«I P = ^^[-sin^fCos^fFoK.m? ) - sin 2 (9 6 cos 2 ^ (m^<) 
+ cos 2 0~ t cos 2 9~ b F (m| , m~ bi ) + cos 2 0- t sin 2 6 b F (m| , m? 2 ) 
+ sin 2 # t ~cos 2 9 b F (ml,ml) + sin 2 %sin 2 B- h F Q {m\, m\ 2 ) ] . (5) 
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Here mi v m~ h {% = 1,2) denote the stop and sbottom masses, whereas 6^,6^ are the mixing 
angles in the stop and in the sbottom sector, see also Sect. 3.1.1. 

The dominant two-loop correction within the SM of 0(aa s ) is given by [12] 

A P 2 -loop = -APl-loopg— (1 + T /3) . (6) 

It screens the one-loop result by approximately 10%. 

The corresponding corrections of 0(aa s ) in the MSSM have been evaluated in Ref. [24]. 
Contrary to the SM case, these corrections can enter with the same sign as the one-loop 
result, therefore enhancing the sensitivity to the squark effects. 



2.2 Electroweak two-loop corrections to Ap: the gauge-less limit 

The Yukawa contributions of 0(a"j) form a set of leading two-loop contributions entering 
the EWPO only via Ap, where «/ = y 2 /(47r), and y/ is the Yukawa coupling of fermion /. 
For the top and bottom quarks the Yukawa couplings read 

V2m t V2m b 

Vt = — : Wi Vb = w , (7) 

Dsmp v cos p 



where v = \Jv\ + v\. In the SM another subset of leading electroweak two-loop corrections 
to Ap is given by the corrections for large Higgs-boson masses of 0(Gf l Mf I M^ v ) [16]. We 
will focus on the Yukawa corrections in the following. 

In order to evaluate the leading Yukawa contributions of 0(a 2 ) the gauge-less limit has 
to be applied. It consists of neglecting the electroweak gauge couplings gi j2 — > and thus 
also = glv 2 /2 ->■ and Mf = {g\ + gl)v 2 /2 0, while keeping the ratio c w = M w /M z 
and the vacuum expectation value v fixed. Accordingly, Y^ z ,w m e Q- (2) need to be evaluated 
at 0(g 2 2 ) m order to obtain a finite contribution of 0(g° 2 ) to Ap in the gauge-less limit. In 
this limit only diagrams with fermions and scalars contribute to Ap, while no gauge bosons 
appear in the loop diagrams. 

At the one-loop level the only non-vanishing contributions to Ap in the gauge-less limit 
of the MSSM are the fermion-loop and sfermion-loop contributions as given in eqs. (3), (5). 
While the Higgs sector of a general two-Higgs-doublet model yields a contribution to Ap in 
the gauge-less limit, the contribution vanishes once the symmetry relations of the MSSM are 
imposed (see the discussion in Sect. 4.1 below). 

At the two-loop level the gauge-less limit results in the desired Yukawa contributions of 
O(oj). For the quarks and squarks of the third generation this yields in particular terms 
of O^mf/v 4 ) and 0{m\ tan 2 /3/v 4 ). It is easy to see that no other contribution to M w and 
sin 2 9 e fi besides the gauge-less limit of Ap yields terms of this order. 

In the SM the two-loop result for Ap in the gauge-less limit was first obtained for the 
special case M h sm = [14], 

Ap^t|M„.„ = 3jL m . (19-2,2) . (8) 
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This result was then extended to the case of arbitrary values of M h sm [15] . The corresponding 
result is given by 



a SM,a? v Q G 2 J x 2 (10-6x 2 + x 4 ) 2 /l-uA 
Ap 2 _ Ioo * p (M^m) = 3^^777, j-3 ^— w\og [~^) (9) 



x 2 (x 2 -4) , fl-w 
H iu log 1 



2 ° V 1 + w 



( 9 X 4 \ , , on X 2 (10 — 6X 2 + X 4 ) , o / ox 

(6 + 6^ 2 - yj log (x 2 ) + 3 l 2(x2 _ 4) ^ log 2 (V) 



+ 25-4x^ + 71 



2 2 +8 — 6X 2 + X 4 — lOtfX 4 + 6ti7X 6 



wx 



2x 2 (x 2 -4) 



2 10 — 6x 2 + x 4 (\ - w 

+ 07i7X Ll 2 

x J — 4 



3(x 2 -l) 2 (x 2 -2) T . 



X 2 



Li 2 (l-x 2 )}, 



where x = M H sM/m t , and w = a/1 — 4/x 2 . The effect of going beyond the approximation 
M h sm = turned out to be numerically very significant. While the numerical value of the 
result in eq. (8) is rather small due to the accidental cancellation of the two terms in the last 
factor of eq. (8), the result is about an order of magnitude larger for values of M h sm in the 
experimentally preferred region. As a consequence, the result for the 0(a 2 ) corrections to Ap 
with arbitrary Higgs-boson mass as given in eq. (9) provides a much better approximation of 
the full electroweak two-loop corrections to the EWPO [7-10] than the limiting case where 
M h sm = 0, eq. (8). As an example, for M h sm = 120 GeV the resulting shifts in Mw and 
sin 2 # eff are -10 MeV and +5 x 10~ 5 , respectively. 

Within the MSSM also the contributions involving the bottom Yukawa coupling can be 
relevant at large tan/3. The corresponding contributions of O(af), 0(a t a b ), and 0(al) to 
Ap have been obtained in Ref. [27] in the limit of heavy scalar quarks. In this limit only the 
top and bottom quarks and the Higgs bosons (and Goldstone bosons) of the MSSM appear 
in the loops. The results turned out to be numerically relevant, leading to shifts in Mw 
and sin 2 9 c s of up to 12 MeV and 6 x 10~ 5 , respectively. Since in the gauge-less limit the 
couplings of the light CP-even Higgs boson of the MSSM to fermions become SM-like, the 
0(a 2 ) correction in the MSSM can be well approximated by the corresponding correction 
in the SM, as given in eq. (9). Potentially larger effects compared to the SM case can be 
expected from the contribution of supersymmetric particles (with not too heavy masses), 
since these corrections do not have a SM counterpart. 



3 The 0(oL 2 t ), 0(oL t OL h ), 0(oL 2 h ) contributions to Ap 

The purpose of the present paper is to perform a complete calculation of the 0(a 2 ), 0(a t ctb), 
and 0(a 2 ) contributions to Ap in the MSSM, including the contributions of supersymmetric 
particles. This means that all diagrams have to be evaluated (applying the gauge-less limit) 
that contain top and bottom quarks, their scalar superpartners stop and sbottom, and Higgs 
bosons or higgsinos. 
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The contributions to Ap at 0(a^), 0(a t at,), 0(af) can be grouped into three classes (see 
Fig- 1): 

(q) diagrams involving t/b quarks and Higgs bosons (see also Ref. [27]), 

(q) diagrams with t/b squarks and Higgs bosons (see Fig. 2 for generic diagrams), 

(if) diagrams with higgsinos (containing also quarks and squarks) (see Fig. 3 for generic 
diagrams) . 

The generic diagrams shown in Figs. 2, 3 have to be evaluated for the Z boson and the 
W boson self-energy. 



Q 
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V V 
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Figure 2: Generic Feynman diagrams of class (q). V denotes either W or Z, q is either a £ 
or a b, and 0, x denote Higgs and Goldstone bosons. 

In the following sections we describe the necessary ingredients for the evaluation of these 
contributions, starting with the relevant sectors of the MSSM and then turning to the tech- 
nical evaluation of the two-loop diagrams and the counterterms. 



3.1 The relevant MSSM sectors 

Here we specify the MSSM contributions that are relevant for the O(af), 0(a t a b ), 0(al) cor- 
rections. As explained above, the calculation involves the gauge- less limit where M w , M z — > 
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Figure 3: Generic Feynman diagrams of class (if). V denotes either W or Z, q is either a £ 
or a b, while q is a £ or a 6, and if denotes a higgsino (neutral or charged). 





V V V V 



Figure 4: Generic Feynman diagrams for the vector boson self-energies with counterterm 
insertion. V denotes either W or Z, q is either a t or a 6, while g is a £ or a 6 

(keeping c w = Mw/Mz fixed). Accordingly, we discuss the implications of the gauge-less 
limit for the different sectors of the MSSM. 
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3.1.1 The scalar top and bottom sector 

We start with the stop and sbottom sector. In order to fix the notation, we explicitly list 
their mass matrices in the basis of the interaction eigenstates £l,£_r and &l, &r: 

w2 _ ( M i +^ 2 + cos2/3(i-|4)Ml m t X t \ 

_ V m tX t Ml+m 2 t + lcos2(3s 2 w M 2 z J ' [W) 

U 2 _ ( Ml+ml + cos 2(3 {-\ + \sl)M z m b X b \ 

M ~ b ~ [ m b X b Ml R +ml-\cos2PslMl J W 

where 

m t Xt = m t (At — n cot (3 ), m b X b = m b (A b — //tan/3 ). (12) 

Here A t denotes the trilinear Higgs-stop coupling, A b denotes the Higgs-sbottom coupling, 
and /j, is the higgsino mass parameter. 

In order to obtain the mass eigenvalues the mass matrix has to be diagonalized. This is 
done with the help of the unitary matrix U^, 

WM\U^ = diag(r4, ml), q = t,b, (13) 

which we will choose real since we consider the MSSM without CP-violation. 
Furthermore, SU(2) gauge invariance requires the relation 

M? l = Ml (14) 

between the soft supersymmetry-breaking parameters in the left-handed components of the 
squark doublet. Consequently, there are only five independent parameters in the t/b sector. 
The masses and mixing angles are connected via the relation 

El^lX = £l^| 2m l +ml-m 2 t -M 2 z clcos2f3. (15) 

i=l,2 i=l,2 

In the gauge-less limit the terms proportional to M§ in the diagonal entries of the mass 
matrices and in eq. (15) vanish. 

Except where stated otherwise, we will assume universality of all three soft supersymmetry- 
breaking parameters in the diagonal entries of the stop/sbottom mass matrices, 

Msusy = M h = M iR = M bR . (16) 

The common squark mass scale is denoted as M SUSY . 



3.1.2 The Higgs sector 

Contrary to the SM, in the MSSM two Higgs doublets are required [31]. At the tree-level, 
the Higgs sector can be described with the help of two independent parameters (besides g± 
and #2 ): the ratio of the two vacuum expectation values, tan/3 = i^/fi, and Ma, the mass 
of the CP-odd A boson. The diagonalization of the bilinear part of the Higgs potential, i.e. 
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the Higgs mass matrices, is performed via orthogonal transformations with the angle a for 
the CP-even part and with the angle (3 for the CP-odd and the charged part. The mixing 
angle a is determined through 



tan 2a = tan 2(3 



Ml ■ Mi 



7T 



M\ 



Ml 



One gets the following Higgs spectrum: 



2 neutral bosons, CP = +1 
1 neutral boson, CV = — 1 
2 charged bosons 
3 unphysical Goldstone bosons 



< a < 



h,H 
A 

H + ,H~ 
G°,G + ,G- 



(17) 



(18) 



At the tree level, the Higgs boson masses expressed through M z , M w and Ma are given by 

1 



Mi = 



Ml 



Ml + Ml- (Ml + M 2 ) 2 - AM 2 A M\ cos 2 2(3 
M 2 A ■ M'-y ■ \J (Ml + M|) 2 — AM 2 A M z cos 2 2(3 



(19) 

(20) 
(21) 



We use the Feynman gauge, which implies the mass parameters M 2 G = Mf , M 2 ± = for 
the unphysical scalars G° and G ± . 



In the gauge-less limit the Higgs sector parameters satisfy the following relations: 



Ml 



sin a = — cos (3 , cos a = sin (3 , 



. 



(22a) 
(22b) 
(22c) 



In the gauge-less limit furthermore the mass of the lightest CP-even Higgs boson vanishes 
at tree level, 

M h = . (23) 

Because of the accidental cancellation in the SM result for M h sm = 0, see eq. (8), it is 
desirable to retain the dependence on M^ as much as possible in the MSSM result. In this 
way the numerical impact of the M/j-dependence can be studied, which within the MSSM is 
formally a higher-order effect. This is particularly interesting in view of the fact that higher- 
order corrections to the masses and mixing angles in the MSSM Higgs sector are sizable, see 
e.g. Ref. [32] for recent reviews. 

We will therefore discuss the implementation of the gauge-less limit in some detail. In 
particular, we will investigate in how far a consistent result for Ap can be obtained without 
imposing eq. (23). We will also briefly discuss the case where eq. (22b) is relaxed, see Sect. 4 
below. For higher-order corrections in the Higgs sector we use the results as implemented 
into the code FeynHiggs [33-36]. 
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3.1.3 Higgsinos 



In the gauge-less limit the contributions from the chargino and neutralino sector reduce to 
those of the higgsinos. 

The diagonalization matrices for the chargino and neutralino sectors in this limit are 
given by (see Ref. [2] for our notation) 



U = V = 




1 



N = 



V2 



( \ 



1-1 

\ o o i i / 



(24) 



and the corresponding elements of the diagonalized mass matrices are 

m~± = (0, +/i), m^o = (0, 0, -fi). 



(25) 



All entries corresponding to gauginos are zero since the gaugino couplings vanish in the 
gauge-less limit. Note that the negative sign in m^o has to be taken into account; the 
physical masses of the charged and neutral higgsinos are all equal to in the gauge-less 
limit. 



3.2 Evaluation of the Feynman diagrams 

The amplitudes of all Feynman diagrams shown in Figs. 2,3,4 have been created with the 
program FeynArtsS [37], making use of the MSSM model file [38]. Dirac algebra and traces 
have been evaluated using the program TwoCalc [39], and the reduction to scalar integrals 
has been performed in two ways using either the package TYReduce [40] or the routines 
built into TwoCalc, based on the reduction method of Ref. [39]. As a result we obtained the 
analytical expression for Ap depending on the one-loop functions A and B [41] and on the 
two- loop function Ti 34 [39,42]. For the further evaluation the analytical expressions for A , 
B and T134 have been inserted. 

In addition to the two-loop diagrams, one-loop counterterms corresponding to the renor- 
malization of divergent one-loop subdiagrams have to be taken into account. The whole 
calculation can be performed both in dimensional regular izat ion [43] as well as in dimen- 
sional reduction [44]. Since no gauge bosons appear in the loops, both regularization schemes 
preserve gauge invariance and supersymmetry for the present calculation. Therefore the nec- 
essary counterterms correspond to multiplicative renormalization of the parameters in the 
MSSM Lagrangian. 



3.3 Counterterms 



The renormalization constants that are relevant in the gauge-less limit are 

5m t , Sm b , <5m|, 5m|, S9 h 5m^, bm\^ 59~ b , 

5M\, <5tan/3, 5t h , H , 5fi, (26) 
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corresponding to the renormalization of the fermion masses, the parameters of the i/b sector, 
the Higgs sector parameters and tadpoles, and the /z-parameter. It is not necessary to 
introduce wave function renormalization constants for the fermions and scalar fields since 
they drop out in the sum over all diagrams. 

There are two possible ways to obtain the counterterm contributions: 

1. Generate and evaluate one-loop diagrams with insertions of counterterm vertices, as 
depicted generically in Fig. 4. In order to generate these diagrams, the required coun- 
terterm Feynman rules had to be added to the FeynArts MSSM model file. In the 
explicit evaluation of the counterterm diagrams it turned out that the renormalization 
constants 8M\, 5tan/3, 8t hjH , S/j, corresponding to the Higgs/higgsino sector, drop 
out. Only the quark and squark mass and mixing renormalization constants contribute. 

2. The renormalization transformation Aj — > A, + 5\i for each parameter A, appear- 
ing in eq. (26) is performed directly in the one-loop result Ap^ M 1 + s J JSY (A i ). The 
counterterm contribution for the two-loop calculation is then obtained by expanding 
A Pi-ioop USY (>i + 5\i) to first order in the SX h where the contributions to Ap™+susY 
have been given in eqs. (3), (5). In this setup it is obvious that the renormalization 
constants 5M\, 5 tan. (3, Sth,H, do not contribute, since the one-loop result in the 
gauge-less limit consists only of the quark and squark loop contributions and there- 
fore does not depend on the Higgs-sector parameters. Accordingly, the counterterm 
contributions for the two-loop calculation, Ap ct , can be written as 

A Pct = Y ( Sm f d ™f + E Sm lH + E S 4 d u?) A Pi M ioo S p USY (27) 
f=t,b V i=l,2 1 ij=l,2 V 

In order to have a non-trivial check of the counterterm contributions, we implemented them 
using both approaches and found agreement in the final result. 

Due to supersymmetry and SU(2) gauge invariance, see eq. (14), there are only five 
independent parameters in the t/b sector, leading to eq. (15). As a consequence, not all 
the parameters appearing in eq. (15) can be renormalized independently. Choosing 
as the dependent parameter, its counterterm Srri-^ can be expressed in terms of the other 
counterterms. In the gauge-less limit the relation reads 

K = [4p(E ra 2 K - i^il 2 K +2«K -O s <2 



i=l,2 



- 2 U\ x u\ x (m? x - ml 2 )5u\ 2 - 2m t 5m t + 2m b Sm^j , (28) 
where the renormalization transformation of the mixing matrix is defined as 

Ufj ^ Uf, + SUf, , SUf^SuiUl, 8u% = -5u%. (29) 

In order to define the renormalization constants one has to choose a renormalization 

scheme. For the SM fermion masses m ti b we always choose the on-shell scheme. This yields 
for the top mass counterterm 

5m t = \m t [Re S ti (m t 2 ) + Re S tj? (m?) + 2Re S ts (m 2 t )] , (30) 
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with the scalar coefficients of the unrenormalized top-quark self-energy, in the Lorentz 

decomposition 

E t (p) = ^_£ ti (p 2 ) +0u + E tR (p 2 ) +m t ^ ts (p 2 ) , (31) 

and analogously for the bottom mass counterterm (in order to take higher-order QCD cor- 
rections into account, we use an effective bottom quark mass value of mj = 3 GeV). For 
the five independent t/b sector parameters we choose either the on-shell [45,46] or the DR 
scheme. The precise definitions will be given in the following section. 

4 Renormalization prescriptions and result for Ap 

As explained above, the strict implementation of the gauge-less limit in the evaluation of 
the 0(af), 0(a t ab), 0(af) contributions to Ap would imply that the mass of the lightest 
CP-even Higgs boson has to be set to zero, see eq. (23). In the SM case, where M H sm is a 
free parameter, it turned out that the two-loop Yukawa contribution to Ap yields a much 
better approximation of the full electroweak two-loop corrections to the EWPO for (realistic) 
non-zero values of M h sm than in the limit M h sm = 0. It is therefore of interest to investigate 
the impact of non-zero values of also for the MSSM, where is a dependent quantity 
that is determined by the other supersymmetric parameters. 

It has been observed already in Ref. [27] that the pure fermion contributions of class (q) 
(see Fig. 1) may consistently be obtained even if eq. (23) is not employed. In this section 
we discuss this issue in detail and explain the physical origin of this behavior. Based on 
this result we show how the calculation of all three classes of contributions to Ap, i.e. (q), 
(q), and (H), can be organized in such a way that Mh can be set to its true MSSM value 
essentially everywhere. We will use the resulting expression in order to study the numerical 
effect of non-zero values for the new corrections calculated in this paper, namely the 
squark contribution (q) and the higgsino contribution (H). 

4.1 Higgs sector 

In order to discuss the implementation of the gauge-less limit it is useful to compare the 
MSSM case with the one of a general two-Higgs-doublet model (2HDM). In the following 
"2HDM" is to be understood as a two-Higgs-doublet model including squarks and higgsinos, 
but without any supersymmetric relations imposed on them. The MSSM can be regarded 
as a special case of a 2HDM, with supersymmetric relations for couplings and masses. In 
the 2HDM without these coupling relations, Ap is also well-defined and can be calculated 
at 0(a 2 ), 0(a t ab), 0(af) in the gauge-less limit. The corresponding two-loop diagrams are 
identical to the diagrams of the classes (q), (q), (H) in the MSSM. However, in contrast to 
the MSSM, the Higgs-boson masses in the 2HDM are independent parameters and do not 
have to obey eqs. (22), (23) in the gauge-less limit. 

The essential difference between the MSSM and the 2HDM case concerns the renormal- 
ization and counterterm contributions. Restricting ourselves in a first step to class (q), these 
contributions can be decomposed in the MSSM and the 2HDM as 

ApSsm = Ap£ oop + Ap£ ct (32) 
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Aa&m = Ap^ioop + ApLlt + ApSL , (33) 

respectively. Here the two-loop contribution Apj,^, and the counterterm contributions 
from the t/b doublet Ap///_ ct are identical in the two models, while the Higgs sector coun- 
terterm contribution Ap^_ ct appears only in the 2HDM result. 

As mentioned above, in the MSSM there are no one-loop contributions from the Higgs 
sector and correspondingly no Higgs sector counterterm contributions at 0(a 2 ), 0(a t ab), 
0(al). In the 2HDM, the Higgs sector one-loop contribution to Ap reads 

a 2 r 

A ^2HDM i/ rr ( aj2 n/T 2\ 

Ap ^ = T28^[ Fo(M ^ Ma) 

+ sin 2 ((3 - a) (F (Ml±, M 2 H ) - F (M 2 A , M 2 H )) 

+ cos 2 03 - a) (F (M 2 H± ,M 2 h ) - F (M 2 A , M 2 )) ] (34) 

in the gauge-less limit. Note that this contribution indeed vanishes if the MSSM gauge-less 
limit relations (22) hold. The counterterm contribution from the Higgs sector at the two-loop 
level can be obtained from this expression as 

Ap ( £l ct = I SM 2 d M 2+5t a npd tan p + 5smad sina \Ap 2U ^ p . (35) 

\<f>=h,H,A°,H± J 

Since F (x,y) and d x F (x,y) vanish in the limit x = y, we find that Ap^_ ct vanishes if 
M H = M H ± = M A , cos(/3 - a) = 0, M h = arbitrary (36) 

or 

M H ± = M A , 5M 2 H± = 5 M\, M h ,M H ,a = arbitrary. (37) 

The relations in eq. (36) are the same as the constraints imposed by the gauge-less limit 
in the MSSM except for the fact that Mh = is not necessary. The observation made in 
Ref. [27] that the class (g) contributions to Ap can be evaluated in the MSSM in a meaningful 
way for non-zero values of M h can be understood from eq. (36). For class (q) the two-loop 
diagrams and the fermion sector counterterms are identical in the MSSM and the 2HDM. If 
the relations in eq. (36) hold, Ap^_ ct = in eq. (33), so that eq. (32) and eq. (33) become 
identical. Thus, the calculations in the MSSM and the 2HDM are the same in this case. 
This means that the result of class (g) derived in the MSSM for non-zero M h is well-defined 
and consistent, as it corresponds to a certain special case of the general 2HDM result. 



4.2 Inclusion of the t/b sector in the on-shell scheme 

For the full set of contributions to Ap, also the sfermion diagrams of class (q) and the 
higgsino diagrams of class (H) have to be taken into account. In the following, as explained 
above, we consider a 2HDM including also stops, sbottoms and higgsinos (although without 
any supersymmetric relations). The analogy of the calculation in the MSSM and the 2HDM 
does no longer hold, since the sfermion sector renormalization differs in the two models. 
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As discussed above, supersymmetry and SU(2) gauge invariance imply that not all pa- 
rameters in the squark sector can be renormalized independently in the MSSM. Choosing 
m~ bi as the dependent mass in the MSSM, its renormalization constant 8m? is given by 
eq. (28), and no independent renormalization condition can be imposed on it. We will refer 
to the expression for 8m? in terms of the other counterterms of the fermion and sfermion 

sector as given in eq. (28) as the "symmetric" renormalization, 8m? L I 

o d 7 oi I symm 

In the on-shell scheme [45], the three other squark masses ??Tf 1 2 02 are defined as pole 
masses, and the mixing angle counterterms can be defined via on-shell mixing self-energies: 



8m\ = Re£ /i (m|) 



8u 
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2(m 2 - - m 2 - ) 



for fi = £1,2,62; 



for f = t,b. 



(38) 
(39) 



In a 2HDM with squarks, on the other hand, an on-shell renormalization can be applied 
for all four squark masses. In this case 8m? ^ is given by 



<Hlos 



(40) 



where Eg is the b\ self-energy. Hence there is a mass shift 



Am- = 8m? I — 8ml I ~~ 

01 bi I symm b\ I OS 



(41) 



at the one-loop level between the mass parameter m?^ as given by the "symmetric" renor- 
malization and the physical pole mass. 

For the class (q, H) the decomposition of Ap in the two models is given by 



^PMSSM 



^A^-loop 



^Ptb-ct 



Ap 



tb—ct, symm 

^2HDM — ^2-loop + ^Ptb-ct + ^Pfb-ct, full OS + ^H-ct 



(42) 
(43) 



Here Ap:f' corresponds to the "symmetric" renormalization of the squark sector in 

' to— ct, symm 

the MSSM described above. Ap£' H ^ „ „ ^„ denotes the contribution from the full on-shell 

' tb—ct, full OS 

renormalization of all squarks. As one can see from eqs. (42) and (43) the MSSM result differs 
from the 2HDM result even for the case where Ap^_ ct = 0. The MSSM result therefore does 
not correspond to a special case of the 2HDM expression. 



4.3 Result for Ap in the on-shell scheme 

The total result for Ap at 0(a 2 , a t a b , a 2 ) in the MSSM is given by the sum of eqs. (32) and 
(42), 

ApM'** = Ap$ SSM + A&Z . (44) 

As discussed above, in Apj^g^ the "symmetric" renormalization in the sfermion sector has 
to be applied, leading to the relations (15), (28) for m?^ and 8m?^. The contribution of class 
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(q,H) can be rewritten by using the mass shift as defined in eq. (41) (see also Ref. [24]), 
leading to the expression 

A _ a (?) , f\M^) » o a A SUSY (Ar\ 

LA(J — LAp MSSM -f- L±P MSSMt full OS + LAm b 1 °m?^ lA Pl-loop > 

where Apj^ full os is given by 

Ap£lL, full os = Ap££ p + Ap®2 + A P |f) fuU QS . (46) 

A°mssm ^ s the result for the fermion-loop contributions as obtained in Ref. [27] (employing 
an on-shell renormalization of the fermion masses and inserting non-zero values for Mh). 

(a fr\ 

A°mssm full os * s ^ ne contribution of the squark and higgsino diagrams obtained by renor- 
malizing all sfermion masses, i.e. including m? , on-shell, while the last term in eq. (45) 
is a symmetry-restoring contribution involving Am? . The one-loop sfermion contribution 
A°i U io^ P nas been defined in eq. (5). 

Comparing eqs. (46) and (43) shows that the "full OS" contribution in eq. (46) is UV- 
finite already for the partial gauge-less limit of eq. (36), according to the discussion of the 
previous two subsections. 

Consistency requires that the mass shift in eq. (41) has to be UV-finite as well. One can 
easily check that this requires to take into account both squark/Higgs and quark/higgsino 
loops. Correspondingly, because of the necessity of this shift only the sum of the (q) and 
(H) contributions to Ap is physically meaningful in the MSSM. 

Moreover, the mass shift Am^ is only finite in the gauge-less limit, i.e. it can only con- 
sistently be evaluated if all the gauge-less limit relations (22a)-(22c) and Mh = 0, eq. (23), 
are used. The last term in eq. (45) can therefore only be obtained in the approximation 
where Mh = 0. 

The expression in eq. (45) represents the main result of this paper. For the first term on 
the right-hand side of eq. (45), Ap^ SSM , we keep the full dependence on Mh- As explained 
above, this is possible because this term is not affected by the renormalization in the sfermion 

sector. For the second term, Ap^'ssL full0 s> we wm keep the M/j-dependence as well and 
compare to the strict gauge-less limit case where M h = 0. The mass shift Am? entering the 
last term in eq. (45) is evaluated for M h = 0. 

The last term in eq. (45) can be expressed using eqs. (42) and (43) as 

A™ 2 f) n An SVSY — A rS&'fy I A^Q'H) I (A7\ 

iim~ bi o m 2^ Pl _ loop - ^P fb _ cttSyui J Mh=0 ^Pf b - ct ,f u \ios\M h =o ■ 
Correspondingly the full result can be rewritten as 

A 0(9.9.*) - An (q) + An^ + An %&) + An^ 

^P ~ ^PMSSM + ^2-loop + ^Ptb-ct + ^Pfb-ct, full os 

I [a JP) I _ AjQ'fy I (AS.) 

~ r [_ ^tS-ct, symmlM h =0 ^tb-ct, full OS \M h =o\ ' V °> 

All contributions in the first line of eq. (48) can be evaluated by keeping the full M h de- 
pendence. For the other parameters of the Higgs sector we impose the gauge-less limit as 
specified in (22a)-(22c). 
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As a result of eq. (37), the gauge- less limit can be relaxed in another way. If the sum 
of all contributions (q, q, H) is considered, the relation 5M^ ± = 8M\ in eq. (37) is valid. 
As a consequence, in the evaluation of the first line of eq. (48) it is not even necessary to 
use the gauge-less limit for sin a and Mh- Instead, sin a and Mh can be set to their true 
MSSM values. We will discuss the case where the gauge-less limit is relaxed also for these 
two parameters below. 



4.4 Renormalization in the DR scheme 

As an alternative to the on-shell scheme in the squark sector, we also consider the DR 
scheme. In this scheme the counterterms of the soft supersymmetry-breaking parameters 
are defined to be pure divergences. The squark mass and mixing angle counterterms receive 
finite contributions corresponding to m t) b in the squark mass matrices (10), (11): 

5m|| fin = (yfSMpH).. foif i = i 1 , 2 ,b 1 ,2; (49) 
Su{ 2 L = ± 2 J / 12 , (50) 

"I fin m 2 _ m 2 V I 

h h 

^-L = *»/L ( %' 2 % ) • < 51 > 

The result for Ap^'^ in the DR scheme follows from eq. (48) by replacing Ap^' H ^ 

to — ct, full OS 

by the corresponding counterterm resulting from eqs. (49)-(51). As a consequence, the 
terms in the second line of eq. (48) vanish. The results in the DR scheme depend on the 
renormalization scale /i DR . 



5 Numerical analysis 

In this section the numerical effect of the electroweak two- loop correction eq. (45), or equiv- 
alently eq. (48), is analyzed, using the formulas in eq. (1) to obtain the corresponding shift 
in My/ and sin 2 9 e g. In addition to the full MSSM correction resulting from Ap^ ,H \ we also 
present the effective change compared to the SM result (where the SM Higgs boson mass has 
been set to Mh). This effective change can be decomposed into the contribution from class 
(q) and from classes (q,H). The contribution from class (q), which was studied in Ref. [27], 
reads 

V 9) (MSSM - SM) = Ap^ SSM - Apl M Cl(M H sM = M h ), (52) 

where Ap™^ has been given in eq. (9). The contribution from classes (q,H) is given by 

ApGft = Aft full os + Am?a m , ApSUSY p , (53) 

where Mh = is used in the second term. Here and in the following we drop the subscript 
"MSSM" for simplicity. 

As SM input parameters we use the values m t = 178.0 GeV,rrib = 3 GeV. The bottom 
quark mass is to be understood as an effective bottom quark mass, taking into account 
higher-order QCD corrections. 



16 



5.1 Impact of relaxing the gauge-less limit for and sin a 



In the first step we study the impact of evaluating Ap^<^) M full0 s ( see ec l- (53)) for the true 
value of the lightest MSSM Higgs-boson mass Mh rather than for Mh = 0. Accordingly, we 
compare the effect on the EWPO resulting from Ap^(M h ) + Ap^(M h ) and Ap^(M h ) + 
Ap^(0). 

We have investigated the numerical effect of keeping the dependence on M h in the squark 
and higgsino contributions for various MSSM scenarios. Fig. 5 shows an example where the 
numerical impact on the prediction of Myj and sin 2 # e g- is quite sizable. The EWPO are 
given as a function of Ma with Msusy = — A.t,b — 400 GeV, p = 800 GeV and tan f3 = 50. 
The effect of keeping a non-vanishing value of Mh in the squark and higgsino contributions 
amounts to about +5 MeV in M w and —3 x 10~ 5 to sin 2 9 e s for all considered Ma values. 
The effects for other MSSM scenarios are typically smaller than for the example shown in 
Fig. 5. Unless otherwise stated, we will always keep the full Mh dependence in the results 
shown below. The difference between the result with and without the Mh dependence can 
be employed for estimating the residual theoretical uncertainties from unknown higher-order 
corrections, see the discussion in Sect. 5.4 below. 



-10 



-12 



■■ 400 GeV, A, = A b = -400 GeV, tanp = 50, |x = 800 GeV 




200 



M h = in (q, H) 



300 400 

M A [GeV] 
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A. CO 



Figure 5: AM W and Asin 2 # eff are shown as a function for Ma for the case where the full 
dependence on the mass of the light CP-even Higgs boson is kept, Ap^(M h ) + Ap^ H \M h ), 
and for the case where the strict gauge-less limit for Mh has been applied in the squark and 
higgsino contributions, Ap^(M h ) + Ap^(O). 



Fig. 6 illustrates the numerical effect of relaxing the gauge-less limit on sin a. As discussed 
at the end of Sect. 4.3, the sum of the contributions of classes (q, q, H) can be evaluated in a 
meaningful way even if sin a and Mh are set to their true values in the MSSM instead of their 
values in the gauge- less limit. Since the corresponding shift in Mh is usually quite small [33] 
we do not analyze the effects arising from different choices for M H and use its gauge-less 
value throughout the paper. The situation is different for the Higgs mixing angle a. Here the 
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Figure 6: AMw and A sin 2 9 e s are shown for the case where the Higgs mixing angle a obeys 
either the full tree-level relation, eq. (17), or is fixed by the gauge-less limit, eq. (22b). 



full tree-level value sin a as given in eq. (17) can significantly deviate from its gauge-less 
value, sin a gl = — cos (5. Fig. 6 shows the results for AMw and A sin 2 6 c q based on sin a fu11 
and sino; gl . The parameters are chosen in such a way as to maximize the influence of sin« fu11 
vs. sino; gl . The value tan/5 = 6 is rather small, and e.g. together with Ma = 150 GeV it 
leads to sina fu11 = -0.31 and sina gl = -0.16. For M SUSY = // = 400 GeV and A tfi = 
—800 this parameter set is in agreement with all experimental constraints from Higgs boson 
searches [47,48] and 6-physics [49]. Fig. 6 shows that even in this scenario the numerical 
effect of relaxing the gauge-less limit on sin a is negligible. We have checked that this holds 
in general. In particular for larger tan (3 and/or Ma the effect is even smaller. Therefore we 
will always set sin a to sin a gl in the following. 



5.2 Dependence on supersymmetric parameters 

In Figs. 7, 8 and 10 we explore the numerical impact of Ap^ q ' q,H ' on Mw and sin 2 6 c r for 
various MSSM parameter choices. The values are chosen such that experimental constraints 
are fulfilled for most parts of the parameter space. Fig. 7 shows a scenario with large tan (3, 
tan/3 = 50, and M SU sy = M A = 300 GeV and fj, = 500 GeV. The results are plotted as 
functions of the stop-mixing parameter X t = A t — /i/ tan/3 (see eq. (10)), and we chose 
Ab = A t . The two-loop contributions AMw and Asin 2 ^ are decomposed into the SM 
result, Ap 2 _ 1 ^ p (M ff sM = M h ), as given in eq. (9) (shown with reversed sign for better 

visibility), Ap^(MSSM - SM) as given in eq. (52), and Ap^' ' as given in eq. (53). For the 
latter contribution both the result with the correct MSSM value for Mh and with Mh = 
is shown. We find that Ap^' H ^> induces shifts in M w and sin 2 9 e g of up to +8 MeV in M w 
and —4 x 10" 5 in sin 2 6^. The corrections are significantly larger than the effective change 
compared to the SM arising from class (q), Ap^(MSSM — SM). The impact of relaxing 



18 



the gauge-less limit on M h in Ap^' H ^ is clearly visible, although not as pronounced as in 
Fig. 5. It should be noted that small mixing in the stop sector (in this scenario values of 
\X t \ < 350 GeV) is disfavored by the LEP Higgs searches [47,48], i.e. the dependence on 
Mh is largest where its value is already experimentally excluded. For small values of \X t \ 
the super symmetric contribution Ap^(MSSM — SM) + Ap^' H ^ is almost as large as the 

SM result, Ap2^{^; p (M H sM = M h ), and largely compensates it. For large values of \X t \ the 
supersymmetric contribution reduces the SM result by about 40%. 
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Figure 7: AMw and A sin 2 e ff are shown as a function of X t in a scenario with large tan (5. 
The two-loop contribution involving squarks and higgsinos, Ap^' H \ is shown for the correct 
MSSM value of Mh and for Mh = 0. For the class (q) the effective change from the SM to 
the MSSM is shown and compared with the pure SM contribution (with the sign reversed 
for better visibility). 



In Fig. 8 we show a similar plot for a parameter scenario with small Higgsino mass, 
p = 200 GeV, and tan (5 = 6, M SU sy = 400 GeV, M A = 300 GeV. The contribution of 
Ap^' H ^ amounts to about 1-2 MeV in M\y and — 1 x 10~ 5 in sin 2 # e fr in this case. The 
fermion loop contribution Ap^(MSSM — SM) is very small here because the small value of 
tan f3 does not lead to an enhancement of in the MSSM with respect to the SM. 

Fig. 9 shows the one-loop results, Ap^j^ , corresponding to the scenarios of Figs. 7, 8. 
Due to the larger value of Msusy and the small value of tan (3 the one-loop contributions for 
the second scenario are relatively small. The region of small \X t \ is again ruled out by LEP 
Higgs searches. The largest effects visible in Fig. 9 are thus experimentally excluded. Com- 
paring the one-loop with the two-loop results, one can see that the two-loop contributions 
from Ap^' H ^ amounts to about 10% of the one-loop supersymmetric contributions. 

A common feature of the two scenarios, visible in Figs. 7, 8, 9, is that both the one- and 
two-loop supersymmetric contributions first decrease for increasing \X t \ until a minimum is 
reached in the vicinity of X t ~ — 2Msusy- For even larger mixing one stop mass becomes 
very small and the supersymmetric contributions increase again. 
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Figure 8: AM W and A sin 2 # cff are shown as a function of X t in a scenario with small \i 
and tan/5. The two-loop contribution involving squarks and higgsinos, Ap^' H \ is shown for 
the correct MSSM value of and for = 0. For the class (q) the effective change from 
the SM to the MSSM is shown and compared with the pure SM contribution (with the sign 
reversed for better visibility). 
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Figure 9: One-loop SUSY contributions to AM^ and Asin 2 ^ e fj are shown as a function of 
X t . The parameters correspond to the two scenarios analyzed in Figs. 7 and 8. 
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5.3 Results in SPS scenarios and renormalization scheme 
dependence 

Fig. 10 shows the results for Ap^' H ^ in the SPS la benchmark scenario [29] for a moderate 
value of tan/3 = 10 and four different combinations for p, and M A , 



(/i/GeV, M A /GeV) = (200, 200), (200, 1000), (500, 500), (500, 1000). 



(54) 



In order to display the dependence on the scale of supersymmetry, we start from the nominal 
values of the MSSM parameters corresponding to the SPS la point [29] (besides \x and M A 
that are chosen as specified in eq. (54)) and vary the parameters M S usy and A t>b using a 
common scale factor; the results are then shown as functions of Msusy- 1 The range of Msusy 
values shown in Fig. 10 has been chosen such that compatibility with Higgs-boson mass [48] 
and 6-physics [49] constraints is ensured for most parts of the parameter space. For small 
values of M S usy the corrections differ by up to 4 MeV depending on the choice of M A and \x. 
Smaller values of M A and \i result in larger corrections to Mw and sin 2 # e ff. In all cases the 
result decreases with increasing Msusy as expected. The corresponding supersymmetric one- 
loop contributions induced by ApfY.^ are shown in Fig. 11 for comparison. The two-loop 

correction from Ap^' H ^ amounts up to 25% of the MSSM one-loop result. 
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Figure 10: The shifts AMy/ and A sin 2 9 e g induced by Ap^' H ^ are shown as a function 
of M SUS y in the SPS la scenario for four combinations of and M A = 200, 500 GeV and 
H = 200, 500, 1000 GeV. 



We now study the renormalization scheme dependence of the one-loop and two-loop 
results for three benchmark SPS scenarios. Besides the "standard" scenario SPS la, we also 
investigate the SPS lb scenario, which is characterized by a larger tan (3 value, tan (3 = 30, 

1 More precisely, for the SPS points the soft supersymmetry-breaking parameters ^ ^ for the left- 
arid right-handed t, b are all slightly different. Msusy is identified with M; . 
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Figure 11: The shifts AM W and Asin 2 ^ induced by the supersymmetric one-loop contri- 
butions are shown as a function of M SU sy in the SPS la scenario for /x = 200, 500, 1000 GeV 
and tan (3 = 10. 



and SPS 5, which involves a relatively light t [29]. Fig. 12 shows the one-loop results for 
the three scenarios, while Figs. 13, 14, 15 display the two-loop results. As above, the results 
are shown as functions of M S usy- We have started from the nominal values of the MSSM 
parameters for the three benchmark points and varied the parameters MgusY, A t>b , /i (for the 
DR results also the scale /i DR ) using a common scale factor. The actual SPS la, SPS lb and 
SPS 5 benchmark points correspond to Msusy = 495.9, 762.5, 535.2 GeV, respectively [29]. 

For a meaningful comparison of the results in the on-shell and the DR renormalization 
schemes, the input parameters in the two schemes have to be physically equivalent, which 
implies that they are numerically different. Since the parameters in the SPS scenarios are 
defined as DR parameters, they can directly be used as input parameters in the DR scheme. 
The corresponding input parameters for the calculation in the on-shell scheme are obtained 
by requiring 

(m| + 5m 2 f ) os = (m 2 ~ + 6m 2 f ) m (55) 

Ji Ji Ji Ji 

for the squark masses and similarly for the mixing angles. 

In the one-loop results Ap^j S Y for the three SPS scenarios shown in Fig. 12 the squark 
sector parameters correspond to the on-shell scheme. The shift in the precision observables 
induced by Ap^f^ is found to be particularly large for the SPS 5 scenario, as a consequence 
of the large splitting between the squark masses in this scenario. 

In Figs. 13-15 we show the one- loop result parametrized in terms of on-shell parameters 
(dotted line) and the two-loop (q, H) results obtained in the DR (full line) and the OS scheme 
(dot-dashed line), in all cases relative to the one-loop result parametrized in terms of the 
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Figure 12: The shifts AM W and Asin 2 ^ eff induced by the supersymmetric one-loop contri- 
bution Ap^ u f ^ p are shown for the three benchmark scenarios SPS la, SPS lb and SPS 5 as a 
function of M SUSY . The parameters of the squark sector correspond to the on-shell scheme. 



DR parameters. Accordingly, the three lines in each plot correspond to 



The pure two-loop correction in the DR scheme is given by the full line, while the two-loop 
correction in the on-shell scheme corresponds to the difference between the dot-dashed and 
the dashed line. 

The numerical impact of the two-loop correction Ap^' H ^ in the scenarios SPS la, lb 
amounts to about 5-6 MeV in M w and —3 x 1CT 5 in sin 2 6 eS for small M SUSY and decreases 
to about 1 MeV in M w (-0.5 x 1(T 5 in sin 2 6 eS ) for larger values of M SUS y- For SPS 5 the 
corrections are slightly smaller. While in the scenarios SPS la, lb the two-loop results in the 
two schemes are very close to each other, a larger deviation is visible in the SPS 5 scenario. 
In the latter scenario the two-loop corrections in the on-shell scheme are less than 1 MeV, 
while in the DR scheme they are more than twice as large. Comparison with the one-loop 
results given in Fig. 12 shows that the two-loop corrections amount to about 10% one-loop 
MSSM contribution. 

The comparison of the renormalization schemes shows that the scheme dependence is 
strongly reduced by going from the one-loop to the two-loop level. At the one-loop level, 
where the scheme difference is entirely due to the different input parameters for the squark 
masses and mixing angles, the difference between the on-shell and the DR scheme is of 
0(1 MeV) in M w . Taking into account the two-loop corrections reduces the difference 
below 0.1 MeV for SPS la,b and about 0.2 MeV for SPS 5. 




(56) 



23 



8 



6 — 



CD 



SPS1a, relative to 1-loop DR 
1 -loop OS 
2-loop DR 
2-loop OS 




CO 



500 
M 



600 700 
SUSY t GeV ] 



1000 



Figure 13: AM W and A sin 2 6^ are shown in the SPS la scenario as a function of M S usy- 
The results for the one-loop contribution expressed in terms of on-shell parameters and for 
the two-loop result Ap^ u |^ p + Ap^' H ^ in the on-shell and the DR scheme are given relative 
to the one- loop result expressed in terms of DR parameters, see eq. (56). 
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Figure 14: AM W and Asin 2 # cff are shown in the SPS lb scenario as a function of M SU sy- 
The results for the one-loop contribution expressed in terms of on-shell parameters and for 
the two-loop result Ap^ U [ S ^ p + Ap^' H ^ in the on-shell and the DR scheme are given relative 
to the one- loop result expressed in terms of DR parameters, see eq. (56). 
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Figure 15: AMyy and Asin 2 # c g- are shown in the SPS 5 scenario as a function of Msusy- 
The results for the one-loop contribution expressed in terms of on-shell parameters and for 
the two-loop result Ap^Yf^ + Ap^' H ^ in the on-shell and the DR scheme are given relative 
to the one- loop result expressed in terms of DR parameters, see eq. (56). 



The size of the two-loop corrections for SPS la,b is found to be much larger than the 
difference between the two schemes at the one-loop level, which is only about 1 MeV for 
these scenarios. This indicates that the difference between the results in two renormaliza- 
tion schemes, if taken as the only measure for estimating the theoretical uncertainties from 
unknown higher-order corrections, may result in a significant underestimate of the actual 
theoretical uncertainty. The SPS 5 scenario, on the other hand, is an example where the 
two-loop corrections turn out to be smaller than the scheme difference at one-loop order. 

Finally we compare the two-loop results for the (q, H) contributions obtained in this 
paper with the two-loop QCD corrections of 0(aa s ) as obtained in Ref. [24]. In Fig. 16 we 
show the results in the on-shell scheme for the three SPS scenarios as a function of M S usy 
(as explained above). For SPS la and lb both corrections are roughly of the same size and 
compensate each other to a large extent. Only for the case of SPS 5 the QCD corrections are 
significantly larger than the two-loop Yukawa corrections. Both the QCD and the Yukawa 
corrections are non-negligible in view of the anticipated future experimental accuracies and 
larger than the current theoretical uncertainties in the SM. 



5.4 Estimate of unknown higher-order corrections 

As discussed above, the theoretical evaluation of the EWPO in the SM is significantly more 
advanced than in the MSSM. In order to obtain an accurate prediction for the EWPO within 
the MSSM it is therefore useful to take all known SM corrections into account. This can be 
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Figure 16: The effect of the two-loop Yukawa corrections from squark and higgsino loops is 
compared with the squark-loop corrections of 0(aa s ). AM W and Asin 2 # cff are shown in 
the three SPS scenarios as a function of M SUSY in the on-shell scheme. 



done by writing the MSSM prediction for the observable 0(0 = M w , sin 2 6 e g, . . .) as 

Omssm = S m + Omssm-sm , (57) 

where 0$m is the prediction in the SM including all known corrections, and Omssm-sm is the 
difference between the MSSM and the SM predictions, evaluated at the level of precision of 
the known MSSM corrections. The expression given in eq. (57) contains higher-order contri- 
butions that are only known for SM particles in the loop but not for their superpartners (e.g. 
two-loop electroweak corrections beyond the leading Yukawa contributions calculated in this 
paper and three-loop corrections of 0(aa 2 s )). In the decoupling limit where all superpartners 
are heavy and the Higgs sector becomes SM-like, the result of eq. (57) obviously yields a 
more precise prediction than a result based on only those corrections which are known in the 
full MSSM. In this case the second term in eq. (57) goes to zero, so that the MSSM result 
approaches the SM result with M h sm = M^. For lower values of the scale of supersymmetry 
the contribution from supersymmetric particles in the the loop may be of comparable size 
as the known SM corrections. In view of the experimental bounds on the masses of the 
supersymmetric particles (and the fact that supersymmetry has to be broken), however, a 
complete cancellation between the SM and supersymmetric contributions is not expected. 
It therefore seems appropriate to apply eq. (57) also in this case. 

Expressing the predictions for the EWPO as in eq. (57) implies that the theoretical uncer- 
tainties from unknown higher-order corrections reduce to those in the SM in the decoupling 
limit. In the SM, based on all higher-order contributions that are currently known, the 
remaining uncertainties in M w [8] and sin 2 9 e g [9] have been estimated to be 

8M™ = 4 MeV, 6 sin 2 6^ = 5 x 1(T 5 . (58) 
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Below the decoupling limit an additional theoretical uncertainty arises from higher-order 
corrections involving supersymmetric particles in the loops. In the following we will estimate 
this additional theoretical uncertainty in the prediction of Mw and sin 2 # e g- depending on the 
supersymmetric parameters. We will provide estimates for the uncertainty for three values 
of the squark mass scale, M S usy = 200, 500, 1000 GeV. A similar approach of estimating 
the remaining uncertainties from unknown higher-order corrections with dependence on the 
supersymmetric parameters has recently been applied to the Higgs sector and implemented 
in the program FeynHiggs2.2, see Ref. [36] for details. 

The remaining uncertainties from unknown higher-order corrections involving supersym- 
metric particles mainly arise from the following sources: 

• Electroweak two-loop corrections beyond the leading Yukawa corrections evaluated in 
this paper: 

We estimate the numerical effect of these corrections by assuming that the ratio of the 
subleading electroweak two-loop corrections to the two-loop Yukawa corrections is the 
same in the SM as in the MSSM. Inserting the known SM corrections [7, 8] we infer 
an estimate of the possible size of the missing supersymmetric electroweak two-loop 
contributions. 

• 0(aa s ) corrections beyond the Ap approximation: 

We estimate the size of these corrections by assuming that the ratio of the contribution 
entering via Ap to the full result is the same as for the known SM result [50]. 

• 0(aa 2 ) corrections: 

We use three different methods for estimating the possible size of these corrections. 
The unknown ratio of the 0(aa 2 ) supersymmetric contributions to the 0(aa s ) super- 
symmetric contributions can be estimated by assuming that it is the same as for the 
corresponding corrections in the SM [13] (estimate (a)) and, using geometric progres- 
sion from lower orders, by assuming that it is the same as the ratio of the 0(aa s ) 
supersymmetric contributions and the 0(a) supersymmetric contributions (estimate 
(&)). As a further indication of the possible size of unknown corrections of 0(aa 2 ) 
we vary the renormalization scale of a s {pP Yl ) entering the 0(aa s ) result according to 
iTit/2 < /i DR < 2m t (estimate (c)). It should be noted that this variation of a s (pP R ) 
corresponds to only a part of the higher-order corrections, so that estimates (a) and 
(6) should be regarded as more conservative. 

• 0(a 2 a s ) corrections: 

Similarly as for the 0(aa 2 ) corrections, we again use three different methods for esti- 
mating these corrections. The unknown ratio of the 0(a 2 a s ) supersymmetric contri- 
butions to the 0(a 2 ) (leading Yukawa) supersymmetric contributions can be estimated 
by assuming that it is the same as for the corresponding corrections in the SM [18] 
(estimate (a)) and by assuming that it is the same as the ratio of the 0(aa s ) super- 
symmetric contributions and the 0(a) supersymmetric contributions (estimate (6)). 
As a further indication of possible corrections of 0(a 2 a s ) we change the value of m t 
in the result for the two-loop supersymmetric Yukawa corrections from the on-shell 
value, m° s , to the running mass m t (m t ), where m t (m t ) = mf s /(l + 4/(37r) a s (m t )) 
(estimate (c)). The latter replacement accounts only for a subset of the unknown 
0(a 2 a s ) corrections. 
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• Electroweak three-loop corrections: 

As an indication of the possible size of these corrections we use the renormalization 
scheme dependence of our result for the supersymmetric two-loop Yukawa corrections, 
see Figs. 13-15. 

We have evaluated the above estimates for the three scenarios SPS la, SPS lb, and 
SPS 5, each for M SUS y = 1000 GeV, 500 GeV, and for M SUS y < 500 GeV 2 (as above we 
have varied Msusy, ^t,b and \x using a common scale factor). The estimated theoretical 
uncertainties for My/ arising from the different classes of unknown higher-order corrections 
are shown in Tab. 1. The result given in each entry corresponds to the largest value obtained 
in the three considered SPS scenarios. The three numbers given for the 0(aa 2 ) and 0(a 2 a s ) 
corrections correspond to the estimates (a), (b) and (c) described above. 



Msusy 


<500 GeV 


500 GeV 


1000 GeV 


0(a 2 ) subleading 


6.0 


2.0 


0.8 


0(aa s ) subleading 


1.8 


0.9 


0.5 


0(aa 2 ) 


3.0, 5.3, 1.5 


1.4, 1.1, 0.7 


0.9, 2.2, 0.5 


0(a 2 a s ) 


1.5, 2.2, 1.4 


0.6, 0.8, 0.4 


0.2, 0.2, 0.2 


0(a 3 ) 


0.3 


0.3 


0.3 



Table 1: Estimated uncertainties for M\y in MeV for different classes of unknown higher- 
order corrections involving supersymmetric particles are given for three values of Msusy- 
The estimates have been obtained using the results for the SPS la, SPS lb, and SPS 5 
scenarios. The three entries for the 0(aa 2 ) and 0(a 2 a s ) corrections correspond to three 
different methods for estimating the uncertainties (see text). 

As expected, the estimated uncertainties associated with the supersymmetric higher- 
order contributions decrease for increasing M SU sy- For the 0(aa 2 ) and 0(a 2 a s ) corrections, 
method (c) that accounts only for a part of the higher-order corrections yields in both cases 
the most optimistic estimate. As discussed earlier, by taking into account the true MSSM- 
value of Mh, certain parts of the electroweak corrections, beyond the leading two- loop Yukawa 
corrections, are included in our result. The difference between Ap^' H \M h ) and Ap^ H \0) 
may be interpreted as an estimate of the size of further, not included higher-order electroweak 
corrections. The numerical analysis in Sects. 5.1 and 5.2 shows that this estimate is typically 
smaller than the estimated total uncertainty in Tab. 1. 

We now combine the values given in Tab. 1 into our total estimate of the remaining 
theoretical uncertainties from unknown higher-order corrections involving supersymmetric 
particles. Adopting the largest of the three values for the 0(aa 2 ) and 0(a 2 a s ) as a conser- 

2 The lowest values considered for M SUSY are 200, 300, 400 GeV for SPSla, SPSlb, SPS5, respectively. 
These are the lowest values shown in Figs. 13, 14, 15. For lower values the parameter points are excluded 
by Higgs mass constraints. 
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vative error estimate and adding the different estimates in quadrature we obtain 



5M W = 8.5 MeV for M SUS y < 500 GeV, 
5M W = 2.7 MeV for M SUS y = 500 GeV, 
SM W = 2.4 MeV for M SUSY = 1000 GeV 



(59) 



An analogous analysis of the remaining higher-order uncertainties can also be carried 
out for sin 2 # e ff. Since parts of the missing higher-order corrections to sin 2 9 c q and Mw are 
related to each other, we employ eq. (1) to infer estimates for sin 2 9 e R from our results for 
M w . This yields 



The full theory uncertainty in the MSSM can be obtained by adding in quadrature the 
SM uncertainties from eq. (58) and the SUSY uncertainties from eqs. (59)-(60). This yields 
5M W = (4.7 - 9.4) MeV and 5sin 2 9 cS = (5.2 - 6.7) x 10~ 5 depending on the SUSY mass 



The estimated uncertainties are smaller than the estimates in Ref. [2] (where an overall 
estimate has been given without analyzing the dependence on the supersymmetric param- 
eters), reflecting the improvement associated with the new corrections calculated in this 
paper. 

The other source of theoretical uncertainties besides the one from unknown higher-order 
corrections is the parametric uncertainty induced by the experimental errors of the input 
parameters. The current experimental error of the top-quark mass [51] induces the following 
parametric uncertainties in Mw and sin 2 9 c q 



Sm^ = 2.9 GeV 5M^ ra,mt = 17.5 MeV, 5 sin 2 9 v c ^ mt = 9.4 x 10" 5 . (61) 



This uncertainty will decrease during the next years as a consequence of a further improve- 
ment of the accuracy on m t at the Tevatron and the LHC. Ultimately it will be reduced by 
more than an order of magnitude at the ILC [52]. The accuracy of the theoretical predic- 
tions for M\y and sin 2 9 c q will then be limited by the uncertainty from unknown higher-order 
corrections (for a discussion of the parametric uncertainties induced by the other SM input 
parameters see Ref. [2]). A further reduction of the uncertainties from higher-order SM-type 
corrections (see eq. (58)) and corrections involving supersymmetric particles (see eqs. (59)- 
(60)) therefore seems to be in order to fully exploit the prospective experimental accuracies 
on M w , sin 2 9 e s and m t reachable at the next generation of colliders [52,53]. 

6 Conclusions 

In this paper we have calculated the two-loop corrections of 0(a^), 0(a t af,), 0(al) to the 
electroweak precision observables Mw and sin 2 # e ff in the MSSM. These are the leading, 
Yukawa-enhanced electroweak two-loop contributions; they enter via Ap and arise from dia- 
grams involving SM quarks, squarks, Higgs bosons and higgsinos. While previously only the 




4.7 x 10~ 5 for Mi 
1.5 x 10" 5 for M 
1.3 x 10~ 5 for M 



susy < 500 GeV, 
susy = 500 GeV, 
SUSY = 1000 GeV. 



(60) 



scale. 
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contribution from the diagrams with quarks and Higgs bosons had been known (correspond- 
ing to the limiting case where all supersymmetric particles are infinitely heavy), we have 
evaluated the complete set of Yukawa corrections including the effects of supersymmetric 
particles. 

We have given a detailed account of the theoretical basis of the calculation, focusing on 
the implications of the parameter relations enforced by supersymmetry. In the gauge-less 
limit that needs to be employed to extract the Yukawa corrections of 0{af), 0(a t a b ), 0(al) 
the lightest MSSM Higgs boson mass Mh vanishes. We have studied in how far the true 
MSSM value for can be taken into account in a consistent way. We have shown that the 
result can be expressed in such a way that the M^-dependence, being formally a sub-leading 
effect, can be kept essentially everywhere and we have compared this result with the case 
where the gauge-less limit is sctrictly imposed. 

In our numerical analysis we have put the main emphasis on the new supersymmetric 
contributions involving squarks and higgsinos. We have analyzed the results of the new con- 
tributions as functions of the squark mass scale Msusy, the stop mixing X t and the higgsino 
and Higgs boson mass parameters fi and Ma- For squark masses of about 300 GeV we find 
corrections of typically +4 MeV in and —2 x 10~ 5 in sin 2 6 C &. In certain parameter 
regions, in particular slightly smaller values of M SU sy or small mixing in the stop sector, 
we find shifts up to +8 MeV in M w and —4 x 10~ 5 in sin 2 # efr . For a wide range of pa- 
rameters, the squark and higgsino two-loop corrections increase the corresponding one-loop 
contributions by about 10%. 

Therefore, the class of diagrams with squarks and higgsinos, which has no SM counter- 
part, gives rise to significant deviations from the SM predictions. This is in contrast with 
the contribution of the diagrams involving quarks and Higgs bosons, which can be well ap- 
proximated by the corresponding SM contribution (setting the SM Higgs-boson mass equal 
to the mass of the lightest CP-even Higgs boson of the MSSM). We have compared our re- 
sult for the two-loop Yukawa correction of 0(a^), 0(a t ab), 0(af) to Mw and sin 2 9 e s with 
the 0(aa s ) correction, which is the only other genuine two-loop contribition to M w and 
sin 2 9 cS known in the full MSSM. We find that the two corrections are of comparable size 
and can largely compensate each other for small values of M SU sy (depending on the other 
supersymmetric parameters). 

We have derived our result in two renormalization schemes, the on-shell scheme and 
the DR scheme for the squark sector parameters. Comparing the two-loop results with the 
one-loop result expressed in terms of the parameters of the two schemes shows a significant 
reduction of the scheme dependence. 

We have shown how the known corrections to the electroweak precision observables in 
the SM and the MSSM can be combined such that the currently most accurate prediction 
in the MSSM is obtained. In the decoupling limit, where all supersymmetric particles are 
heavy, the theoretical uncertainty from unknown higher-order corrections reduces to the 
uncertainty of the SM contribution. For non-vanishing contributions of the supersymmetric 
particles an additional theoretical uncertainty arises from unknown higher-order corrections 
involving supersymmetric particles. 

We have estimated the current uncertainty from unknown higher-order corrections involv- 
ing supersymmetric particles for different values of the squark mass scale M SUSY . This has 
been done using geometric progression from lower orders, employing known results for corre- 
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sponding SM corrections, investigating the renormalization scheme dependence, varying the 
renormalization scale, and taking into account formally subleading M/j-dependent contribu- 
tions. For a squark mass scale below 500 GeV we obtain an estimated uncertainty of about 
8.5 MeV in Mw and 4.5 x 10~ 5 in sin 2 9 e g . These uncertainties reduce to about 2.5 MeV 
in M w and 1.5 x 10 -5 in sin 2 ^ e g for Msusy = 1 TeV. They can be combined quadratically 
with the theory uncertainty from unknown higher-order SM contributions to obtain the full 
MSSM theory uncertainties. While currently these uncertainties (for Msusy < 500 GeV) are 
about a factor of two smaller than the parametric theoretical uncertainties induced by the 
experimental error of the top-quark mass, their impact will become more pronounced with 
the expected improvement of the experimental precision of m t . The new two-loop corrections 
evaluated in this paper have been important to reduce the theoretical uncertainties to the 
present level. Further efforts on higher-order corrections in the MSSM will be necessary in 
order to reduce the theoretical uncertainties from unknown higher order corrections within 
the MSSM to the level that has been reached for the SM. 
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